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Trace maps of two-letter substitution rules are investigated with spe- 
cial emphasis on the underlying algebraic structure and on the ex- 
istence of invariants. We illustrate the results with the general- 
ized Fibonacci chains and show that the well-known Fricke character 
I(x, y, z) = x 2 + y 2 + z 2 — 2xyz — 1 is not the only type of invariant 
that can occur. We discuss several physical applications to electronic 
spectra including the gap-labeling theorem, to kicked two-level systems, 
and to the classical ID Ising model with non-commuting transfer ma- 
trices. 



1 Introduction 

Trace maps have been introduced to calculate the spectrum of ID Kronig-Penncy 
and related models on non-periodic structures like the Fibonacci chain ^]. In 
the sequel, these maps have found a broad variety of applications, ranging from 
electronic and phononic structure to ID scattering and transport (cf. || [|, [|, 
and references therein), from a spin in a magnetic field to kicked two- level systems 
J?], ||] , and from classical to quantum spin systems |], [it , 11 . Within only a couple of 



years, the literature on these subjects thus became widespread and hard to survey. 

On the other hand, there is a considerable amount of activity on the mathemat- 
ical aspects of trace maps, ranging from the systematic treatment of substitution 
rules |l2| , through computer algebraic investigations p3| to powerful analytic and 
algebraic results 14, [15], |l6l , |l 7, 18 pj| , and to the treatment of trace maps as 



dynamical systems as such jj^Of 2l|. Unfortunately, several if not many of the ex- 
act results have not found their way to the explicit treatment of physical systems, 
partially because they escaped notice and partially because they require some non- 
trivial mathematics like, e.g., the general gap labeling theorem for Schrodinger 
operators jlf], 22 1. 



It is this gap between the theory of trace maps and their application which we 
would like to bridge, and we hope that the present article might prove useful for 
that. It is important to find a systematic, but simple formulation of the trace maps. 
By this we mean a formulation that is close to that of the Fibonacci case (being 
the best known) but brings no artificial complications for the generalizations (like 
singularities induced by coordinate systems). It turns out that the formulation 
of jl^] is very appropriate. Here, all two-letter substitution rules (and we will 
only consider those) lead to trace maps that are mappings F from C 3 to C 3 where 
the component mappings are polynomials with integral coefficients only. No poles 
are present as in alternative formulations with rational functions (cf. Jll|), and 
invariants can be described more easily. Therefore, we will follow this path and 
occasionally link it to other formulations. 

Let us now briefly outline how the remainder of this article is organized. In 
Sec. 2, we start with some basic definitions and properties of trace maps, using the 
algebraic approach via the free group of two generators |TJ, [l7j . We then formulate 
the invariance and transformation properties of the Fricke character I(x,y,z) — 
x 2 + y 2 + z 2 — 2xyz — 1 , a proof of which is given in the Appendix. We illustrate 
the results with several examples and discuss, in Sec. 3, the generalized Fibonacci 
chains in some more detail. 

Sec. 4 deals with the application to electronic structure, where we apply the gap 
labeling theorem to the generalized Fibonacci class. For a subclass of models, we 
also present an intuitive, direct approach which does not rely on the abstract result 
and, additionally, gives an indexing scheme of the gaps of the approximants. In Sec. 
5, we very briefly describe the application to a kicked two-level system which follows 
a quasiperiodic SU(2) dynamics. Sec. 6 then deals with the classical ID Ising model 
with coupling constants and magnetic fields varying according to the non-periodic 
chain. Here, the trace map gives rise to a direct iteration of the partition function 
which allows at least an effective numerical treatment of the case of non-commuting 
transfer matrices. This is followed by the concluding remarks in Sec. 7 and an 
Appendix where we prov e the general transformation law of the Fricke character 
I{x,y,z) as given in Eq. (2.16). 



2 Trace maps and their properties 

To explain what a trace map of a two-letter substitution rule is, and what its main 
properties are in simple but general terms, it is the best to start from the free group 
generated by two generators a, b, i.e., from 

F 2 :=«a,b». (2.1) 

Its elements consists of all possible finite words w = w(a,b,a~ 1 ,b^ 1 ). F 2 is an 
infinite group, multiplication is formal composition of words, and the empty word, 
e, is the neutral element. 

A substitution rule g now is a mapping of F 2 into itself, where - for obvious 
reasons - we are only interested in homomorphisms, i.e., in q's with the property 



g(w 1 w 2 ) = g(w 1 )g(w 2 ). 



(2.2) 



Here, the multiplication is the group multiplication in i 7 ^- To specify a homomor- 
phism, it is therefore sufficient to give the images of the two generators, 



g : 



It',, 
Wb 



(2.3) 



where w a and Wb are words in a, 6, a _1 ,6 _1 , for details we refer to [24| |T7j] . The 
concatenation of two homomorphisms, 02 ° 0i, is well defined. However, one has 
to insert the rule of Q\ into that of 02 which is a bit against the usual practice. 
Therefore, we follow jul and define the multiplication 



0102 = 01-02 



02 01- 



(2.4) 



With this multiplication, the set of homomorphisms becomes a monoid with the 
identity as neutral element. We call this monoid 2 = Hom(F 2 ). 

A special role is played by the invertible homomorphisms 0, i.e., by automor- 
phisms. They form a group $2, that is also finitely (but not freely) generated 

n n n 

$2:=«U,a,P» (2.5) 
where the generators are defined by 



w 


U(w) 


o~(w) 


P(w) 


a 


ab 


a- 1 


b 


b 


b 


b 


a 



For the relations amongst the generators, we refer to p. 164 of p4j . 

Automorphisms transform the group commutator K(a, b) :— aba~ 1 b~ 1 to a con- 
jugate either of itself or of its inverse 



g(K(a,b))=K(g(a),g(b)) 



[K(a,b)] 



±1 



(2.6) 



where w is an element of F2 and {K(a, b))^ 1 = K(b, a). This can directly be verified 
by means of the generators, but has also been shown independently p7f . As we will 
see, this does not extend to g $ $2- 

For the characterization of a substitution rule, it is useful to define the corre- 
sponding substitution matrix. First, one maps the elements of F2 to lattice points 
in the square lattice, Z 2 , by simply adding the powers of a and the powers of b 
separately in the word considered. This powercounting induces a homomorphism 
from $2 = Aut( J F 2 ) to G1(2,Z) = Aut(Z 2 ) @ and, even more, from the monoid 
9 2 = Hom(F 2 ) to Hom(Z 2 ) = Mat(2,Z), the latter being the monoid of 2x2- 
matrices with integral elements. Note that Gl(2, Z) consists of precisely those ma- 
trices R S Hom(Z 2 ) with det(i?) = ±1, and it is this determinant which determines 
the exponent in Eq. (2 



We call R e the substitution matrix of the substitution rule 0, if we can read the 
number of a's and b's in the words w a and Wb rowwise, i.e., 



R„ 



#a{w a ) 
#a{w b ) 



#b{w b ) 



(2.7) 



Note that often the transpose of R g is used. With our definition of the product ga, 
see Eq. (p^3|), we have indeed 

Rga = Rg ' R<71 (2-8) 

and, e.g., the Fibonacci rule reads a — > b 1 b — > ba and leads to the matrix 

Rg = ( J I ) ■ ( 2 -9) 

Of course, different substitution rules can result in the same substitution matrix 
and the kernel of the homomorphism from 82 to Mat(2, Z) contains all inner auto- 
morphisms of $2: if g acts as g(a) = waw^ 1 and g(b) — wbw~ l , one has R e = 1. 
The kernel of the mapping from $2 to Gl(2, Z) consists of precisely these inner 
automorphisms |25| while there are further elements in the general case. 

So far, we have only set up the algebraic background. To come to the trace maps, 
we have to interpret the substitution in the free group as a substitution — and thus 
a recursion — in the group S1(2,C) of 2x2 complex unimodular matrices. We simply 
insert two matrices A, B for the letters a, b and continue. But S1(2,C) matrices have 
very special properties as a consequence of the Cayley-Hamilton theorem of linear 
algebra, namely 

A 2 = tr(A)A-l 
A + A- 1 = tr(A)l (2.10) 

A n = Un-^A-Un^x)! 

where x = ^tr(A) and U n (x) are Chebyshev's polynomial of the second kind p8| 
defined by: 

U-!(x) = 0, U (x) = l 

U n+X {x) = 2xU„(x)-U n - 1 (x) . y - > 

For later use, we run the recursion also backwards and formally obtain polynomials 



U n for all n 6 Z. Also, Eq. (2.1C) is then valid for all integers n. 
We now define 

x=\tx{A), y=^tv(B), z=±ix(AB). (2.12) 

Then, as a result of Eq. ( |2 . 1 0| ) and of the general identities tv(A+B) = tr(A)+tr(B) 
and ti(AB) = tr(BA), we find that the traces after one iteration of an arbitrary 
substitution rule g can be expressed by x, y, z again [^2|. This is the corresponding 
trace map 

x \ / f e (x,y,z) 

F e : ( V ) - ( 9 e (x,y,z) | (2.13) 
h g (x,y,z) 

where fg,g g ,h g £ Z[a;,y, z] are three polynomials with integral coefficients. All 
trace maps defined that way are thus C°° mappings ofC 3 into itself - and interesting 
dynamical systems as such. Furthermore we find 



± go ± g w ± a , 



(2.14) 



compare also |14j . The mapping from the substitution rules to the trace maps is a 
homomorphism between monoids. For the mapping from $2 to Diff(C 3 ), the kernel 
is known, compare |2~1| . The mapping induces also a homomorphism from G1(2,Z) 
to the trace maps, the kernel of which is {±]l}- 

One of the most celebrated properties of trace maps is the invariance of the 
polynomial 

I(x, y, z) = x 2 + y 2 + z 2 - 2xyz - 1 (2-15) 
under all trace maps that stem from an automorphism (l], [l3|, |l4|, O. In our 



formulation, it follows either immediately from the Nielsen theorem (Eq. (2.6)) 
(observing that A G S1(2,C) implies tr(A) = tr(A -1 )) or, alternatively, from the 
generators U, a, P in Eq. ( |2.3| ) and their trace maps by simple substitution into 
Eq. fl2~T5p : 



U 



P 



Ra 



1 1 

1 



-1 
1 



1 

1 



F„ 








This remarkable property appeared in the Fibonacci case in Q , but the theory of 
Fricke characters, where I(x,y, z) belongs to, is much older (cf. [jl7| and references 
therein): It goes back to the last century j23|. 

But this is not the end of the story because one can show that the Fricke charac- 
ter follows a specific transformation law under any substitution rule, not just under 
automorphisms. Indeed, for g G Hom(F2), one finds 



I(F g (x, y, z)) = P g (x, y, z) ■ I(x, y, z) 



(2.16) 



where P e (x,y,z) € Z[x,y, z] is a polynomial with integral coefficients, called the 
transformation polynomial of g. This relation was conjectured on the base of com- 
puter algebraic manipulations and proved shortly after In the Appendix, 
we give an independent and slightly more complete proof of it. 

This relation, together with the general properties mentioned above, has a num- 
ber of strong consequences. The reader is invited to find some himself. We list the 
following: 



• The variety A4q = {(x, y, z) €C 3 | I(x, y, z) = 0} is an invariant surface for 
all trace maps. 



• (1, 1, 1) is a fixed point and {(1, — 1, —1), (—1, 1, —1), (—1, —1, 1)} is an invari- 
ant set of all trace maps. 



• If (0, 0, 0) is not a fixed point of F e , its image lies on Mo- 
rn If q G $2, (0,0,0) is a fixed point of F B . Conversely, F e (0,0,0) = (0,0,0) 

implies P e (Q, 0, 0) = 1 and Vf e (0, 0, 0) = (0, 0, 0), but not P e = 1. 

• P e (0,0,0) is either or 1. 

• If P e is constant, we have either P e = 1 or P Q = 0. 

• If g is of finite order, g€ $2 and P g = 1. Furthermore, one has g n = e with 
n = 1, 2, 3, 4, or 6 and F™ = Id with m = 1, 2, or 3. 

• The substitution g is invertible if and only if P g = 1 . 

• The substitution g> has nontrivial kernel if and only if P e = 0. 

• The substitution g is injective but not onto if and only if P e ^ const. 

• The set {F B \ g € $2} of invertible trace maps is a group isomorphic to 
PG1(2,Z). 

We cannot give the proofs here, many of which can be found in the work of 
Peyriere and coworkers, compare p9[ and references therein. Let us only remark 
that one takes profit from using trace orbits of SU(2) matrices and some standard 
arguments from calculus. 

Instead, we would like to remark that, for all trace maps, the dynamics is com- 
pletely integrable on A4q. Any triple (x, y, z) can be obtained as traces of diagonal 
matrices, cf. Appendix. Those in turn can be written as diag(exp (ia),exp (— ia)), 
a £ C, wherefore the multiplication is easily done by counting the exponents of the 
matrices involved. For any approximant, this can be expressed by the eigenvalues 
of the substitution matrix R g and is a straightforward generalization of the cor- 
responding property of the Fibonacci chain. This way, one sets up an interesting 
relation to the theory of pseudo Anosov maps, see 

H Hi for details. 

3 The generalized Fibonacci chains and new in- 
variants 

In this chapter, we will illustrate some of the properties in the class of generalized 
Fibonacci chains. The latter are defined by the substitution rules 

where k,£ G Z. If the context is clear, we will frequently drop the index (k,£). R e 
is unimodular if and only if k = ±1. This coincides with g being an automorphism, 
i.e., invertible. From det(H e ) = —k we also see that k ^ gives non-singular cases. 



What are the corresponding trace maps? A one-page calculation results in 

F e : I y\ I g(x,y,z) J (3.2) 

\ z J \ h{x,y,z) J 

with g(x, y, z) = U k -i{x)Ut-\{y) ■ z - U k -i(x)Ui- 2 (y) ■ x - U k - 2 {x)Ui-x{y) ■ y + 
Uk-2{x)Ut- 2 {y) and h(x, y, z) = g(x, y, z)\i^t+i- 

This is equivalent to the formulation of pi] . There, the iteration is given as 
a difference equation of the form y n +\ = r(y n ,y n —i,y n — 2 ) with rational function 
r. The latter has singularities as a pure consequence of the coordinates chosen. 



Therefore, Eq. (3.2) is advantageous. Another page of calculation establishes the 



transformation polynomial to be 

P e (x,y,z) = (E4-i(:r)) 2 . (3.3) 

The cases k = 1 and k — —1 yield P g = 1 (and hence the invariance of I(x,y, z)) 
which we know already from the invertibility of g in these cases (independently 
of £). The substitution rule g has nontrivial kernel for k = (e.g., the kernel 
contains a~ e b) wherefore P e = is the necessary consequence. The remaining 
cases (k £ {— 1,0,1}, £ £ Z) belong to the class of injective, but non-invertible 
substitutions. Many properties of Sec. 2 can therefore be demonstrated within this 
class of examples. 

One interesting question is the existence of further invariants. That this indeed 
is possible can be seen from the case I — 1, k = 2 where we find 

y \ 

2xz-y . (3.4) 

4xyz - 2x 2 - 2y 2 + 1 J 

As can be checked directly, this leaves the polynomial H invariant, 

H(x,y,z) = (4x 2 - l)y-2xz (3.5) 











• 


- ( 









This could be expected from the close relationship of Eq. (3.4) to the period doubling 



map pi, 13, 32| where an analogous invariant exists. H is different from / and does 
in fact foliate the single invariant sheet {I(x,y,z) = 0}, see Fig. 1. 

A closer look shows that the ID lines on {/ = 0} are an immediate consequence 
of R g having an eigenvalue of modulus one. The motion of the dynamical system on 
{I = 0} is integrable and related to torus maps. But if there is a direction neither 
contracted nor expanded, one could expect another invariant of the trace map F e 
which foliates at least the surface {/ = 0} and, eventually, even the whole C 3 . 

Now, our substitution matrix R e has eigenvalues A± = \ {I ± \J i 2 + Ak) . If they 
are rational numbers, they are automatically integers. This happens if and only if 

k = m-e + m 2 (3.6) 



with m integer, where we find A+ = £ + m and A_ = — m. Eq. (3.4) corresponds 



to m = 1, and it turns out that it is not a singular case. Indeed, for k = £ + 1 in 



Eq. (3.2), we find the amazingly simple invariant 



H(x, y, z) = U e+1 (x) - y-Ui(x) ■ z. 



(3.7) 



This is a generalization of Eq. ( |3.5| ), for a proof one needs the identities Uf - 
Ui-iUi+i = 1 and UfUi-i — Ue+iUg-2 = U±. Also for the case m = —(£ + 1) 



Eq. (3.7) gives the right answer, only A+ and A_ are interchanged. 



The remaining cases are m = — 1 and m = 1 — 1. Both result in fc = 1 — £ and are 



equivalent. From the extension of the recursion, Eq. (3.2), to negative indices, we 
know that f7_( n+2 ) = — U n for n G Z. By means of the identity Uf_ 1 — U1U1-2 = 1 
one can show the invariance of the polynomial 

H(x, y, z) = U e _i(x) ■ y - U^ 2 (x) ■ z = -(U^ x (x) ■ y - U^(x) ■ z). (3.8) 



This again is a very simple expression and resembles that of Eq. (3/7). We think 
that no further invariants exist in the class of models investigated here but could not 
finally settle this question. We can combine the two examples shown formulating 
k = 1 ± £ and H^'(x, y, z) = ±([/ ± ^ +1 )(ir) • y — U±t{x) ■ z). This covers all cases 
in the class of generalized Fibonacci chains where the corresponding torus map has 
a unimodular eigenvalue. 

Let us close this section with a short remark why invariants are important. 
The mappings are three dimensional and could show the dynamics of 3D discrete 
dynamical systems. If an invariant exists, they locally look like 2D systems and 
inherit the well understood properties of them like period doubling routes to chaos 
etc., compare the Appendix of pj| . 

Now, the Fricke character is an invariant precisely for the trace maps of au- 
tomorphisms, but many interesting trace maps do not stem from automorphisms. 
Nevertheless, period doubling can occur (as in the so-called period doubling map). 
The existence of invariants like H or H can explain why — and a further investi- 
gation is in progress. 



4 Electronic spectra and gap labeling 

Amongst the applications of trace maps to physical systems, that to electronic spec- 
tra is perhaps the most widespread one ||, |3). The transfer matrices attached to 
the different intervals describe the propagation of the amplitudes of the wave func- 
tion through the chain. In the case of the continuous Schrodinger equation, the 
transfer matrices for plane waves belong to SU(1,1) || |3q| . The corresponding 
trace map acts in IR 3 and allows the description of the spectrum, while the under- 
lying matrix system can be used to identify the states. The best way to proceed is 
to tackle the full Schrodinger equation 

h 2 

( A + V(x))V(x) = E^(x) (4.1) 

2m 

where V{x) is the potential on the chain. But things are a little bit complicated 
because several quantities become energy dependent. Especially the very invariant 
I(x, y, z) shows a significant dependence on the energy of the incident wave. There- 



fore, most authors restrict themselves to the tight-binding case, where Eq. (4.1) is 
replaced by the discrete version (in suitable units) 



(4.2) 



where I(x, y, z) — \{V2 — V\) 2 . V\, V% are the potentials on the two different inter- 
vals. Although one misses some typical features of the continuous case ]5|j, many 
generic properties of the spectra remain the same, e.g., in both cases many spectra 
are Cantor sets with zero Lebesgue measure (see [ fl6| for a criterion). The price 
one has to pay for going from the full Schrodinger equation to the tight-binding 
approximation is that in the latter case one needs a whole one parameter set of 
Hamiltonians which correspond to the single Hamiltonian in the Schrodinger equa- 
tion @. 

Nevertheless, looking at the gap labeling of the spectra, we can make the re- 
striction to the tight-binding case for simplicity without loss of information. We 
emphasize that all results can be extended to the continuous equation, compare 



For the tight-binding Hamiltonian of Eq. (4.2), the transfer matrices read: 



and, therefore, 

x = itr(Ti) = \{E-Vx) 

V = £tr(T 2 ) = \{E-V 2 ) (4.4) 

z = itrCZVTO = \[E-Vx){E-V2)-\. 

One could now think of a specific substitution rule g and investigate the spectra by 
iterative techniques. 

Here, we want to follow a slighty different path. First we want to present the 
famous gap labeling theorem formulated by Bellissard and coworkers. Because it 
uses non-trivial mathematics we will interpret it in simple words and apply the more 
concrete version to the well-known case of the generalized Fibonacci sequences. 
Afterwards, we will give an intuitive way of understanding how the gap labeling 
theorem works. 

Bellissard and coworkers proved the following theorem using techniques of C* 
algebras and K-Theory 



Let H be a Hamiltonian as in (4.2) with a bounded potential and A^ the 
C* algebra generated by the family of operators obtained from H. by translation. 
The possible values of the Integrated Density of States (IDOS) on the gaps of the 
spectrum are given by the image of the Kq group of the C* algebra Ay^ of the 
Hamiltonian Tt under the trace per unit volume, t*(Ko(A?jj)). 

This is a general result no matter what the dimension D is. Thinking of the 
ID case (fh^) with a potential that takes on only finitely many different values, 
the gap labeling theorem tells us that the possible values of the IDOS on the gaps 
are given by all possible frequencies of all possible words in the infinite chain of 
the substitution g. But one can express all these frequencies of words with length 
greater than one by the frequencies of the words of length one and two. This leads 
to the "concrete" gap labeling theorem mM: 



Let 7i be a Hamiltonian as in (4.2) with the potential given by a primitive 



substitution, g, on a finite alphabet, A. The possible values of the IDOS on the 



gaps in the spectrum are given by the Z(A _1 ) module (resp. the part contained in 
[0,1]), constructed by the components of the normalized eigenvectors and in 2 ) 
to the maximal common eigenvalue A of the substitution matrix for one-letter words 
and for two- letter words Mf\ 
The substitution g is called irreducible if, for any pair of letters a, 6 from the 
alphabet A, the word g k (a) contains 6 for some k. If k can be chosen independently 
of the letters a, 6, then g is called primitive. To guarantee the existence of a (half- 
infinite word as a fixed point of g, one considers a substitution which generates 
an infinite chain from every letter of the alphabet A. There must be at least one 
letter, b € A say, so that g(b) begins with b, and this letter b must appear in every 
possible chain. 

Let us now apply this theorem. For the calculation of and we need 

the substitution of all two- letter words. Let g(w)— a aia2---a n be a substitution 
of a word w, then (with the usual notation, cf. jl6|) 

Qn(w) = {a ai...a N -i)(a 1 a2...a^)...(a m -ia m ...am+N-2) (4.5) 

is the substitution of an TV-letter word (m being the total length of g(ao) obtained 
by the powercounting described in Sec. 2). In particular, the two-letter substitution 
is: g2(w) = {aoai){a\ci2)---{a-m-ia-m) ■ Our example is the recursion 



a — > o 
b -> b e a k 



(4.6) 



which leads to the matrix sytem: T n +i — • T„, compare |5[ |Tl). The matrix 

M^p is nothing but the transpose of the substitution matrix R s , i.e., 



M (D = ( k ) m ^ = l {A1) 

e I i £ J ' mfc = k + £ 



with eigenvalues 



A± = ~(i± V^ 2 + 4fc). (4.8) 

Let us introduce D = k(k + £ — 1) and exclude the case D = (g is reducible for 
k = and does not increase the length of any word for k + £ = 1). Then, we can 
write the normalized eigenvector to A + as 



Here, normalization is a statistical one, i.e., JwWi = 1. Also, is easy to 

calculate. We find 



aa — » 66 g(aa) — » (66) 

a6 -> 6 m a fc g(a6) -> (66) 

6a -> 6V6 _ ' g(6a) -> (66) f - 1 (6a)(aa) fc - 1 (a6) 

66 -> b e a k b e a k g(bb) -► (66) £ ~ 1 (6a) (aa) fe ~ ^afe) 



(4.10) 



which means 



M (2) = 

Q 



/ ft- 1 
1 
1 

V i i e-i 



ft-1 \ 
1 
1 

£-1 / 



The eigenvalues are 



A e {0,0,i(£± V^ 2 + 4fc)}. 
The (statistically) normalized eigenvector to A + reads 

/ (k + £)(k-l 



(4.11) 



(4.12) 



„(2) = 



+ 
+ 

V -(2ft + £) 



- (*-l)A+ \ 
A + 

+ (fc + l)A+ J 



(4.13) 



because of 1/A+ = (A + — i)/k, and X\ = k 



Now the calculation of the frequency module has to be done. This requires 
some care. First, we observe that the components of can be obtained by inte- 
gral linear combinations of components of v^ 2 \ and the latter are all of the form 
jj[(m((fc + £) - A+) + n(-(2k + £) + (ft + 1)A + )], m,n G Z. This can be rewritten 
as 

— (/I + PA+), (4.14) 
but now with the additional constraints 



p, + v = 0(D) 
p+(k + £)P = 0(D) 



(4.15) 



This is necessary to guarantee m, n € Z in the preceeding expression. Remember 
that 1/A+ = (A+ — £)/k. It is easy to check that multiplication with (A + — £) leads 



to new numbers p! , v' which also fulfil Eq. (4.15). Consequently, the Z(A + X ) module 
turns out to be 



its , a \ \ r 1 A + ^A+ 



/i + v =0 (D) 
fi+{k + i)D = 0(D) 



}. (4.16) 



for the continuous Schrodinger equation and r* (Kq (Af^ ) ) n [0,1] for the tight- 
binding case. 

One has to keep in mind that these are only all possible values the IDOS can 
take. The gap labeling theorem does not tell us whether all gaps are really open. 
The most prominent example where gaps are closed systematically is the Thuc- 
Morse sequence [57|. But there, the eigenvalues A± are integers, and one has to 
check whether the frequency module really requires all linear combinations of the 
components of and or not. The analogous question arises for the generalized 



Fibonacci chains with Eq. (3.6) as condition, but we cannot go into details here. 

The special subclass of the metallic Fibonacci sequences (k = 1,1 € IN) has the 
module 

{-(fx + i>X + ) | ftiieZ, ii + D = 0(£)}, (4.17) 



whereas the largest eigenvalues of the substitution matrices are given by the metallic 
means: 

X+ = ^(t+y/F+l) = \t;l, £,£,...}. (4.18) 

Looking at these chains one is able to give a more heuristic but intuitive way of 
understanding how the gap labeling theorem works. 

Let us introduce the generalized Fibonacci numbers, fn, by 

/Si=W£ ) +/£i, ^=0, f? ] = I- (4.19) 

This choice of initial conditions will prove useful shortly. We observe that two 
consecutive numbers are coprime because the greatest common divisor obeys 
gcd(tf\f n %) = gcd(fi e \£ ■ fi e) + fi 1 !,) = gcdtf^,/^), i.e., coprimality is 
inherited from the initial conditions. 

If we now consider the series of periodic approximants starting from 6, i.e., 
b — > b e a — > (b e a) e b etc., we know that these approximants are optimal in the follow- 
ing sense: a given approximant of length gn \ where 

(4-20) 

is closer to the limit structure than any given other approximant of length L < g n . 
One therefore expects that the structure of the IDOS of this approximant is as close 
to the limit IDOS as possible with approximants up to that length. In particular, 
the values of the IDOS on the gaps converge rapidly. 

It is tempting to select series of IDOS plateaus in the sequence of approximants in 
such a way that the corresponding gaps belong to each other (by quantum numbers, 
strength, symmetry, etc.). Then, if one can find a labeling, the limit n — > oo should 
reproduce the result of the general theorem described above. The chains with the 
metallic means facilitate this procedure. Here, gaps in consecutive approximants 
can be attached to one another by the following simple rule. Given a gap in the 
n-th approximant. Then, in the step from n to n + 1, one always grabs the closest 
IDOS value possible. If this is not unique (which happens only occasionally) one 
takes both possibilities as two branches. One branch will actually correspond to a 
new series of gaps, but that is hard to see in general and does not matter for the 
limit. 

Now comes the trick: the gaps in the n-th approximant lead to IDOS values of 
the form 

1W ^ W , 0<m</W + /W lt (4.21) 

Jn + Jn-l 

which is an exact result from Bloch theory for the periodic approximants. But, m 
can be written as 

m = fi-fW+vfW 1 , li^eZ, (4.22) 
because f n and are coprime. Taking {[i,v) as label, it turns out that 



Eq. (4.22) selects series of IDOS plateaus — and hence gaps — that follow the 



rule mentioned before! Furthermore, it is easy to see that neither possible gaps are 



missed nor index pairs are missing (as long as \i ■ /„ + v ■ € [0, fi^ + f„li\)- 
It now is straightforward to calculate the limit points which gives the set 



f A+-l , , 1 N 



M^eZ}n[o,i]. 



(4.23) 



A similar argument was recently given for the original Fibonacci ch ain 



this procedure is most simple. A simple subs tituti on shows that Eq. ( 4.25 ) precisely 



581 where 



reproduces the result obtained above in Eq. ( 4.17 ), including the modulo condition. 
The labels obtained from Eq. ( 4.23 ) are slightly handier as can be seen in Fig. 2, 
where, as an example, we show the IDOS of the first approximants of the 'octonacci' 
chain, (k — 1,£ — 2), and the index pairs starting with the approximants of 

length 3 and 7. 

Let us finally remark that the coincidence of the abstract and the concrete 
approach has a strong consequence: As far as we calculated, all gaps of the ap- 
proximants involved here are open (for suitable choices of potentials). But from the 
second approach, it is meaningful to conjecture that all gaps are open for the class of 
metallic Fibonacci chains. For a proof, one could try pertubative arguments along 
the lines of |3l| because these chains are quasiperiodic with embedding dimension 
2, compare H. This will be given elsewhere. 



5 Some properties of kicked two-level systems 

As already mentioned in the Introduction, matrices and trace maps derived from 
two-letter substitution rules find applications in a variety of problems. Here, we 
briefly describe the case of SU(2) dynamics. The latter appears naturally in the 
solution of Schrodinger's equation for an electron or another particle with spin 
s = 1/2 in a magnetic field B . 

ifi^ = n B Bat> (5.1) 

with Bohr's magneton fiB — Pauli's spin matrices <Tj, a spinor ^ of two 

components, and 

£> — ( B 3 B 1 -iB 2 \ 

Ba= (B 1+ iB 2 -B 3 l ^ 



We have written the time-dependent equation in Eq. (5.1) because the stationary 
case with constant magnetic field is completely integrable [^] and hence less inter- 
esting than the truly time-dependent one. Here we are interested in a special class 
of time-dependent fields that allow the application of the recursive transfer matrix 
techniques, where tran sfer is now in time. 

Rewriting Eq. (|]l]) slightly, we obtain the structural core as 

* = -m(t)y (5.3) 

with hermitian and traceless operator H(t). The (formal) solution is 



*(t) = U(t,to)9(to) 



(5.4) 



with ^ 

U{t,t Q ) = T[exp(-i / H(t)<1t)}. (5.5) 

The time evolution operator U is a time-dependent SU(2) matrix, and T denotes 
time-odering fl40| . 

Let us now think of a time sequence following the generalized Fibonacci chain, 
cf. (Ul). With U n := U(t n ,to) we obtain 

U n+1 = U k n _ x ■ U l n (5.6) 

for the SU(2) dynamics. The treatment of this matrix system gives full information 
but goes beyond the scope of the present article. For a summary of what can be done 
with that we refer to fl7fl . Here, we follow and consider the trace systems only 
which gives insight into the behaviour of certain observables. Let us formulate it for 
the Fibonacci case [k = I = 1), the extension to the other examples is immediate. 

To be explicit in the formulas, we will consider <5-kicks on the time intervals, 
more precisely, at the beginning of them. Any SU(2) matix A can be written in 
the form 

A = exp (—iancr) = cos (a)l — i sin (a)nd (5-7) 
with a S IR + and n a unit vector (hence, (ncr) 2 = II). From the relation 

[na,n'cr]=2iekimnkn'((T m (5.8) 

we can conclude that transfer matrices of the elementary intervals certainly com- 
mute if the kicks are in the same direction. But then, the real traces run on {/ = 0} 
and stay bounded, \x\ < l,\y\ < 1, \z\ < 1. The orbit belongs to the pseudo Anosov 
system, compare pQ , The question now is what happens off this surface. 

Consider the kicks 5(t)cijrij(T, j = 0, 1 for the two elementary time intervals. 
We find 

— exp (— i / 5(t)ajrij(rdt) — cos (a,j)l — i sin (aj)rij(T (5-9) 
Jo 

with 

(riocr) ■ (riier) = (riorii)l + i(rio x ni)a. (5.10) 
We can work out the product — U^'U^ 1 ' and obtain 

= {cos (a ) cos (ai) — sin (ao) sin (ai)noni}H (5-11) 
+i{no x ni sin (ao) sin (ai) — rig sin (ao) cos (ai) — ni cos (ao) sin (ai)}<7 

This gives the initial conditions for the traces as 

zo = cos(a ), Zi=cos(ai), (5 12) 

^2 = cos (ao) cos (ai) — sin (ao) sin (ai)(noni). 

With some further calculations one can establish the formula 



I(x ,xi,x 2 ) = [(n ni) 2 - 1] • (sin(ao)sin(ai)) 2 



(5.13) 



which shows — 1 < I < 0, as it must be for SU(2). 

One can now start from I = — 1 and go to I = by an adiabatic change of 
parameters. This way, one could try to follow (for the Fibonacci case, say,) the 
period doubling route to chaos, described recently in It is an interesting ques- 
tion whether one could experimentally find the conservative feigenvalues predicted 
theoretically. 



6 Classical ID Ising model with non-commuting 
transfer matrices 

Let us consider a linear chain of N Ising spins <jj £ {±1}, j — 1,...,N, with the 
energy of a configuration a = (o"i, cr 2 , . . . , ct/v) being given by |4]{] 

N 

E (v) = -J2(Jj,3+iVj°'j+i + Hja-j), (6.1) 

3=1 

where we choose periodic boundary conditions, i.e., <Jn+i — The canonical 
partition function of this system is obtained as a sum over all possible config- 
urations a as follows 

N 

Z (N) = J2l[cxp (K^+^a^+hjaj) (6.2) 

a j=l 

with Kj,j + i = Jjj+i/ksT and hj = Hj/ksT. Here, T is the temperature and ks 
denotes Boltzmann's constant. The corresponding free energy per site F^ N ^ is given 
by F( N ) = — j* logZ^ where we absorbed the factor l/(fcgT) into the definition 
in order to obtain a dimensionless quantity. 



Eq. (6.2) can be rewritten in the following way (see |41| ) 

a 

= tr (t^tW-.-TW) (6.3) 

where the elementary transfer matrix ' is given by the 2 x 2 matrix |I) 

T ,„ +1) = ( e*p(K iJ+1 + *±^) exp(-^, +1 ) \ 

\exp(-i^ +1 ) exp(x jJ+1 -^±i) ) ' 

while denotes the transfer matrix of N sites. Note that the transfer matrices 

at different sites in general do not commute and therefore cannot be diagonalized 
simultaneously 

So far we have not specified the coupling constants Kj j + i and the magnetic 
fields hj. A particularly interesting case is that of coupling constants Kjj + \ = 1 



and quasiperiodic fields hj . The ground state properties have recently been analyzed 
exactly p2| . Here, we now want to consider non-periodic chains where Kjj + i and 
hj are modulated according to the two-letter substitution rule g^ k ^ ( |3.l| ) but taking 
only two possible values. To be more precise, consider a particular word w n (wq = a, 
wi = b, w n+ i = Q {k ^(w„ )) wi th length N n = /ft* 5 + fc/ftf, where the /ft^ are 



defined (generalizing Eq. ( 4. IS )) by 



/ft? = + */£? If = 0, = 1 ■ (6-5) 

To u)„ we associate an Ising spin chain with N n = N„ ' sites in the following way: 
We choose K jjj+ i = K® and hj = (i = 0, 1) where i = if the j th letter 
in the word w n is an a and i = 1 otherwise. Hence the chains we consider are 
characterized by four real parameters and few, i = 0, 1, and the two numbers 
k, £ £ Z which determine the substitution rule. Denoting the transfer matrix of the 
chain associated with the word w n by T n , one has the following recursive equation 
for the transfer matrices 

T n+1 = T l n T k n _ x (6.6) 

with initial conditions 

_ / exp^W+fcW) exp(-^W) \ 

1 ~ V cx p ex P \ K{i) - h(l) ) ) 

corresponding to the words w = a and W\ = b. 



Now, we are almost in the position to use the trace map (3.2) to obtain a 
recursive relation for the partition function and hence for the free energy of our 
chains. In order to do so, we define unimodular matrices T '„ = -^-T n where 
d n = Vdet T n and use the trace map ( |3.2| ) for 

x n = itr(T B ) = JLz„ (6.8) 

which yields a recursion relation for the partition function Z n = Xv{T n ) of the 
chain associated with the word w n , compare p3] for an equivalent formulation with 



a nested iteration. The determinants are simply given by (see Eq. (6.5)) 

d n +i = d n d n -i 

= (do^idi)^ (6.9) 

with di — y / 2smlT(2i?Wy, i = 0, 1. In particular, all d„ are non-zero if <io and d\ 
are different from zero, i.e. if both coupling constants i = 0,1, do not vanish. 



Thus the trick we used to be able to apply the trace map (3.2) to the transfer 
matrices of the Ising chain was to split the recursion relation (6.6) for the transfer 
matrices into two parts: one for the trace of the transfer matrices and a rather 
trivial one for their determinant. 



Explicitly, the trace map (3.2) becomes 




Vn I — > ( Vn+i ) = [ g(x n ,y n ,z n ) I , (6.10) 

hyX n ^ Urn Zji) 



compare Eq. ( |3.2| ), with starting values 

tr(T ) _ expif^cosh/iW 
2d o (2sinh(2^(°))) 1/2 



i'0 



tr(Ti) expK^coshh^ 

yo = ~~ — = ~, 7771/2 ( 611 ) 

2d i (2sinh(2X( 1 ))) 1 



ZQ = 



tr(ToTi) 
2do(ii 

exp(ir<°> + if (1) ) cosh(feW + feW) + cxp(-(A:( ) + if W)) 
2 (sinh(2if(°)) sinh(2if «)) 1/2 

Iterating the trace map thus corresponds to performing the thermodynamic limit 
n — > oo. 



In this way the trace map (3.2) provides an efficient way to compute the free 
energy in the infinite size limit numerically. Unfortunately, an analytic solution 
to the recursion relation is only known for the case that the starting point lies on 
the invariant variety Mq = {(x,y, z) e C 3 | I(x,y, z) = 0} which turns out to be 
equivalent to the requirement that the transfer matrices Tq and T% commute. For 
this case, the solution is trivial since one only counts the relative frequencies of the 
two different couplings in the chain. 

It should be added that recursion relations for the magnetization (one-point 
function) can be derived by similar means. So far, we have not been able to extend 
this procedure to higher correlation functions. A renormalization group analysis 
(see |Q) of the field- free case (h^ — /jW — 0) gives the substitution matrix as 
linearization at the zero-temperature fixed point and hence the thermal exponent 
a = 1 as one would expect. Further work along these directions is in progress. 



7 Concluding remarks 

Having described structure and applications of trace maps of two- letter substitution 
rules, we will briefly address further developments. It seems that the mathematical 
aspects arc understood fairly well. Nevertheless, the trace maps as dynamical sys- 
tems on their own right show interesting features that deserve further exploration: 
orbit structure, nontrivial symmetries, reversibility, period-doubling, and pseudo 
Anosov structure, compare for some results. 

The physical applications could still profit from simply soaking up known ex- 
act results. A prominent example is the gap labeling where the combination of 
algebraic methods, generalized Bloch theory and pertubation theory improves the 
understanding of spectra and wave functions. For the latter, the extension to ma- 
trix maps should also be studied in more detail. A similar remark holds true of the 
other applications. Kicked two-level systems and spin systems are not yet brought 
to the understanding possible. 

Another, rather restrictive point is the treatment of two-letter rules: more than 
that is desirable. Several results are known for n-letter substitution rules, compare 



Jl(], H |i| Q , but the generality of the findings is lacking. Many simple prop- 
erties do not extend, e.g., trace maps are higher dimensional and investigation of 
invariants is much more complicated. Only the gap labeling theorem still applies in 
full generality, wherefore further work in this direction will be needed. 
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Appendix 

If F g is the trace map of a two- letter substitution rule g E Hon^i 7 ^), then there is 
a uniquely defined transformation polynomial P g E Z[x,y, z] such that the Frickc 
character I(x, y, z) = x 2 + y 2 + z 2 — 2xyz — 1 shows the transformation law 

I{F e (x, y, z)) = P g (x, y, z) ■ I(x, y, z) (A.l) 

Here, the standard basis is used, i.e., x = ^tr(A),y = ^tr(B),z = ^tr(AB) with 
A,B E Sl(2,£). Instead of the complex numbers, one can also use any other 
algebraically complete field. 

For the proof of this remarkable equation, the variety Mo = {(x,y,z) E 
C 3 | I(x, y, z) = 0} plays an important role. We first observe that 

tr(K(A,B)-l)=A-I(x,y,z) (A.2) 

for unimodular 2x2-matrices, where K(A,B) = ABA X B '~ is the group com- 
mutator. This follows straightforwardly from the Cayley-Hamilton theorem. We 
observe next that any triple (x,y, z) E Mo can be seen as traces of diagonal ma- 
trices since we can simply choose A — diag(a, a -1 ), B — diag(/3 e , f3~ £ ), with 
a = x + \/ x 2 + 1, j3 = y + \Jy 2 + 1, and e = ±1. Then, we find from the condition 
I{x,y,z) = (a quadratic equation in z), the relation 2z — a/3 £ + l/(a/3 £ ), i.e., 
z = ^tr(AB). The two choices of e cover the two possible solutions of the quadratic 
equation. 

Observing that [A,B] = implies K(A,B) = 11, we see from Eq. flA.2j ) that 
I(x, y, z) = implies I(F g (x, y, z)) = 0. On the other hand, we know that 

_ I(F e (x, y, z)) _ tr(K(g(A), {B)) - 1) 
I[x,y,z) tr(iv (A, B) - II) 



Thus, P g is the quotient of two polynomials with integral coefficients. The denomi- 
nator, I(x, y, z), is irreducible in Z[x, y, z], and the numerator polynomial vanishes 



at least on the set A4q. From this, we can conclude that the denominator divides 
the numerator. Since this can been seen in an elementary way, we will briefly give 
the argument. 

Consider a polynomial Q ££[x,y,z] that vanishes at least on Afo- If Q is of 
less than second order in x, y, or z, we can only have Q = 0. To show this, assume, 
without loss of generality, Q linear in z, i.e., Q = £(x, y)+rj(x, y)-z with polynomials 
£ and rj. For arbitrary x, y € C we can find z± = Z±{x, y) so that (cc, y, z+) G A4o 
and (x, y, Z-) S A^o- Consequently, Q must also vanish on these points, which gives 

z+-v(x,y) = Z--rj(x,y). (A.4) 

But we can choose (x,y) such that z + (x, y) 7^ z-(x,y) wherefore r)(x,y) = 0. 
Furthermore, there is a whole e-d isc i nC 2 around (x, y) with z + (x, y) ^ z_(x, y). On 



this disc, rj vanishes due to Eq. (A.4). Thus 77 = onC and, from Q(x,y, z + ) = 0, 
also £ = 0. Hence, Q = under the conditions stated. 

Let us now consider Q = I(F g (x, y, z)), which vanishes at least on M.q. If Q is 
only linear (or less) in z, say, the above argument applies and we have Q = — 
giving P e = as well. But, if Q is at least quadratic in z, we can apply Eulcr's 
division algorithm because the ring Z[x,?/, z] is factorial (4?j and the denominator 
starts with z 2 . We find 

Q = P e -I + r (A.5) 

with r a polynomial of degree < 1 in z. But r obviously vanishes at least on A^o, 
wherefore we must have r = according to the above given argument. 

We have thus rigorously and completely established that I(x, y, z) devides 
I(F B (x,y, z)). Now, it remains to be stated that Z[ar, y, z] is facto rial wherefore 
P e is again a polynomial with integral coefficients. This proves Eq. ( |A.l| ). 
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Figures 




Figure 1: The invariant sheet S = {I(x,y,z) = 0} and an orbit of the trace map 
(k = 2, £ = 1) lying also on the second invariant H (x, y, z) which foliates S. 




Figure 2: The approximants of the octonacci chain, k = 1,£ = 2, with length a) 3, b) 
7, c) 17, d) 41 and the gap labels, (/x, v), starting with the 3- and the 7-approximant. 



